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Abstract 

We propose a resolution of the puzzle posed by the discrepancy between the 
value of the pion-nucleon sigma term inferred from pion-nucleon scattering, and 
that deduced from baryon mass splittings using the Zweig rule. We show that 
there is a significant hypercharge-dependent dynamical contribution to baryon 
masses, not hitherto included in the analysis, which may be estimated using the 
scale Ward identity, and computed by solution of the Schwinger-Dyson equation 
for the quark self-energy. We find that the discrepancy is completely resolved 
without any need for Zweig rule violation. 
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1. The Sigma Term Puzzle 

The sigma term puzzle [Q is one of several instances in which experimental data seem 
to disagree with our understanding of the structure of the nucleon based on naive quark 
model intuition. The pion-nucleon sigma term is defined as the nucleon matrix element of 
the light quark mass term in the QCD Hamiltonian: 



a = fh(N\ [uu + dd) \N), (1.1) 

where, ignoring isospin breaking, m = \(m u + mj. It is physically interesting because 
it relates the pattern of chiral symmetry breaking in QCD to the quark content of the 
nucleon 

The value of a can be determined experimentally by using current algebra to connect 
the matrix element (^T|) to the value of the isospin even pion-nucleon scattering amplitude 
at zero momentum transfer |J. Taking several subtleties in the extraction of the matrix 
element from the scattering data into account [|J leads to the value a ~ 45 MeV, with an 
overall experimental and theoretical uncertainty of perhaps up to 10 MeV. 

On the other hand, the sigma term may also be related to the mass splittings in the 
baryon octet |J. If we assume that there are no strange quarks in the nucleon (according 
to the Zweig rule), or more specifically that (iV|ss|./V) = 0, then the sigma term (|1 . 1|) is 
equal to the octet sigma term 

a 8 = m{N\ (uu + dd-2ss) \N), (1.2) 



which in turn is proportional to the nucleon matrix element of the octet portion 7i 8 of the 
quark mass term in the QCD Hamiltonian: 

a 8 = - ^-—(N\H 8 \N) (1.3) 

1 — m s /m 

where 

H 8 = § (m - m s ) [uu + dd - 2ss) . (1.4) 

However, in the quark model the mass splittings are assumed, using first order pertur- 
bation theory, to be due to the SU(3) octet component of the (effective) strong Hamiltonian 
8 ; this leads to the relation 

M B = M + M X (B\Y\B) + M 2 (B\ (I (I + 1) - \Y 2 ) \B), (1.5) 

where / is the isospin operator, and Y the hypercharge. Eq.( |1.5| ), which expresses the 
masses of all the octet baryons in terms of two parameters, is in excellent agreement with 
the data. The matrix element of the symmetry breaking operator 8 between nucleon 
states can then be determined by SU(3) algebra; 

(N\0 8 \N) = Mi - \M 2 

= M A -M E = | M N - | (M H + M s ) . ^'^ 
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Now, if Og is identified with the octet component of the QCD Hamiltonian, namely Tig, 
using the value of the ratio m s /m = 25 ± 5 determined through current algebra |7| 
eq.( |1.3|) determines ag = 25 MeV, with an uncertainty of about 5 MeV. The large per- 
centage discrepancy between the values of a and as is known as the sigma term puzzle. 

It has now become rather fashionable to see this discrepancy as a consequence of the 
failure of the assumption upon which the identification of a and ag is based, namely the 
Zweig rule; one then infers that the strange matrix element (iV|ss|iV) is relatively large. 
More than 300 MeV of the nucleon mass would then be due to the strange quark conden- 
sate, and kaons could condense out in nuclear matter at low densities [Q. Alternatively, 
the discrepancy could be interpreted as a breakdown of first-order perturbation theory, 
i.e., of the linear dependence of the mass splittings on the symmetry breaking operator 
which is used to derive eq.( |1.5| ) and relate the mass splittings to the matrix element ( |1.3| ). 

Both of these alternatives are rather unpalatable, since they are extremely hard to 
reconcile with the successfulness of the quark model mass formulae which are derived using 
the same assumptions. For example, the Gell-Mann-Okubo mass formula (|1.6|) is accurate 
to a few percent. Also, Weinberg has determined the current quark masses assuming 
again a linear dependence of baryon masses on quark masses; assuming that the operator 
ipijj is proportional to the quark plus antiquark number, the mass splittings are given by the 
difference in quark content of the various hadrons. This is essentially equivalent to eq. (|1.5| ) 
with M2 = 0, and is phenomenologically accurate to about 20% |10| [JTJ in agreement with 
the fact that fitting eq.( |1.5|) to the baryon octet spectrum leads to M 2 ~ 

Here we will show that there is a third possibility which successfully resolves the puzzle; 
namely that whereas the mass splittings are indeed given by a linear mass formula according 
to eq.([L5|), the identification of matrix elements of 0% with those of TCs is incorrect due 
to the fact that most of the baryon mass, and thus a substantial proportion of the mass 
splittings, arises dynamically from the trace anomaly. Once this contribution is taken into 
account, it follows that the octet operator responsible for mass splittings differs from Tig 
due to the presence of an additional term, generated dynamically, which may be viewed 
as a nonsinglet gluonic contribution to the mass splittings. Thus the puzzle is resolved, 
but the quark model results, which follow from the assumption of linear dependence on an 
octet operator, are preserved. 

In section II we show that the conformal anomaly equation leads to an exact relation 
between the sigma term and mass splittings, which implies that the identification of the 
nonsinglet matrix elements of 0% and TL& is indeed spoiled by the presence of a dynamical 
contribution to the masses. We then use a Ward identity to argue that this contribution is 
nonsinglet and large enough to account for the observed discrepancy. In section III we test 
this explanation by attempting a computation of the nonsinglet part of the dynamically 
generated mass, solving the Schwinger-Dyson equation for the quark self-energy in various 
approximations. We find that this gives a value for the sigma term perfectly consistent 
with the data, despite large uncertainties in the infrared dynamics responsible for the mass 
generation. Conclusions are drawn in section IV. 
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2. The Trace Anomaly and Dynamical Mass Generation 

In order to discuss the relationship of the mass splittings to the matrix elements of 
Tig we exploit an exact relation between the matrix elements of the sigma term and the 
nucleon masses. Classically the sigma term ( |1 . 1| ) is equal to the divergence of the Noether 
current for scale transformations, which in turn equals the trace of the energy-momentum 
tensor. In the quantized theory, this implies a Ward identity that allows us to relate the 
matrix elements of the trace of the energy-momentum tensor to the masses of physical 
states. 



2.1. The Scale Ward Identity 

Consider the dilation current |12| jp = x v T^ v . This is the Noether current for scale 



transformations; its divergence is equal to the trace of the energy- momentum tensor T^ v . 
In the quantized theory, the trace of the energy momentum tensor satisfies on-shell the 
operator equation [OX 



d»f D = T% = (1 + Tm)^™^ + ^G^G% , (2.1 



where the sum runs over all quark flavors, G£ u is the gluon field strength, (3 = d< ^{^ is 

the beta-function for the strong coupling a s = |-, and 7 m (/i) = ~—j^r^- is the mass 
anomalous dimension. All operators appearing in eq.( ]2.1| ) are renormalized and normal- 
ordered. The last term on the right hand side of eq.( [2.1| ) is due to the conformal anomaly; 
the anomalous dimension r ) m is present because the anomaly term and the mass term in 
eq.( |2.1|) are not separately scale invariant (i.e., they mix upon renormalization) , while the 
energy momentum tensor is (up to surface terms) [p~3[] . Taking matrix elements of eq. ( |2.1| ) 
generates the scale Ward identities of QCD. 

On the other hand, the forward matrix element of the trace of T^ v between hadron 
states is just the mass of the hadron, due to Lorentz invariance and the absence of a 



massless scalar Goldstone boson [[T2| . It follows that the baryon matrix element of eq. (|2.1 



is 

(B\ (1 + lm) EifmWi + GJJfl) = M Bt (2.2) 

which is the desired Ward identity. 

Eq.( |2.2| ) shows that the first-order perturbative expression which equates mass split- 
tings to nonsinglet matrix elements of the mass term is protected by the scale Ward 
identity, and exact up to quantum corrections. These appear in eq.(|2.2j) in two distinct 
instances, namely, the nonzero value of the anomalous dimension 7 m , which effectively 
rescales the current masses, and the presence of the conformal anomaly contribution. We 
can thus separate Mb into a "current" contribution Mg and a dynamical contribution 
Mg; M B = M% + M§ where 

Mb = (5|I>i^<|fl>, M^ = (B\^G^G% + -y^rm^B). (2.3) 
Notice that this is a scale invariant separation, since both M% and M§ are renormalization 



group invariant ||13||. Accordingly, we can separate the mass splittings into current and 
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dynamical contributions by denning M-p, M% and M-j , M® in analogy with (|1.5| ), with 
Mi = Mi + Aff and M 2 = M§ + M 2 D . The same algebra that led to (|L6|) also gives 
(JV|ft 8 |iV) = Mf - ±M 2 C , so that 

AM = (N\0 8 \N) - (N\H 8 \N) = M® - ±M 2 D . (2.4) 

In the conventional argument || , Mg is neglected with the result that AM vanishes, 
and the sigma term is given by eqns ( |1.3| ) and (|1.6[) . This is presumably done on the 
grounds that the term in ( f2.^ ) proportional to the mass anomalous dimension y m is small 
enough to be ignored, whereas that containing the isosing let gluon operator ^^G^G^ 
can be assumed to make no contribution to mass splittings. We will now show that this 
naive assumption is incorrect; away from the chiral limit M® is actually rather large, 
and with opposite sign to M-p, so that the magnitude of (N\TCg\N), and thus of a§, is 
significantly increased. a 



2.2. Estimation of the Dynamical Mass Splitting 

To this purpose, it is convenient to rewrite the matrix element on the left hand side of 
eq. ( |2.2|) in terms of one-particle irreducible physical couplings, analogously to what is done 
in the pseudoscalar case in order to derive the isosinglet Goldberger-Treiman relation [p~5j| . 
This can be done by defining a connected generating functional, W(Sf?, Sa) where 
is the source for the dilation current jp, and Sa are sources for the set of fields $,4 = 

(B, B, (pi = ipiipi, Q = ^^ G^G" ) (which includes the baryons, the scalar quark 
condensates, and the gluon condensate respectively). Also, as in ref.flT5| we define the 
Zumino effective action T(S®, $^() by Legendre transformation of W(Sj? , Sa) with respect 
to the fields $a (but not jp), so that T(S®, $^4) generates diagrams which are one-particle 
irreducible with respect to <& A l - 

Rewriting the Ward identity (|2.2|) in terms of the generating functionals W and V it 



can then be shown HTlH that 



x ^ s 3 w S 3 W 

Mb = (1 + 7m) \ m *5 Sl 5S B 5S B + 5S Q 5S B 5S B (2 ' 5) 



1 Of course a small nonzero value of is present due to the term proportional to the 

mass anomalous dimension. This is actually as large as 30% of Mf , since at the nucleon scale 
7m(l GeV) = 0.27 (to two loop order) p^| . This contribution has (obviously) the same sign 
as Mi and accordingly if taken into account (while neglecting anything else) would make the 
disagreement between the value of a and that of as worse, because in eq.( |1.3| ) a% would be replaced 



by its rescaled value o~g = as- However including such a contribution while neglecting that 

of the gluonic operator is a rather dubious procedure, since the separation of M§ into a "quark" 
and "gluon" piece is necessarily dependent on the renormalization scale. 
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where Abb (A bb ) is the baryon (antibaryon) propagator. It follows that the one particle 
reducible contributions to the matrix elements of the gluon and quark condensates in 
eq. ( |2.5|) must cancel against each other in order to yield the one particle irreducible coupling 
of eq.( |2.6D , namely 

(B\^G^G^\B)°^ + (1 + 7m ) (B^m^B)^ = 0, (2.7) 
or, in the notation of eq. ( |2.3|) , 

(M^) opr + (Mf ) opr = 0. (2.8) 

This shows immediately that it is quite unreasonable to assume M B to be isosinglet given 
that (just as in the pseudoscalar sector, where similar results hold [|15),[[16|]) one would ex- 



pect the one particle reducible contributions to M B to have a significant flavor-nonsinglet 
component. 

More specifically, the one particle reducible matrix elements are presumably domi- 
nated by diagrams where the various operators couple directly to a meson state with the 
appropriate quantum numbers, which then in turn couples irreducibly to the baryon; 

(Mg ) opr = (B\J2 i m i ^ i \B) opi « 2<O|Eimi^i|0a)^-<0a|SB), (2.9) 



where \(j) a ) are the scalar meson states, and the sum over states implicitly includes an 
integration over momenta. We can roughly estimate the order of magnitude of this correc- 
tion by assuming that it is dominated by the exchange of two ideally mixed scalar mesons 
4>i and <f> s , which are respectively an isosinglet pure u and d state, and a pure s state. 
Introducing meson-to- vacuum coupling constants g<j, and meson-baryon couplings g^BB-i 
normalized according to 

(0|f (uu + dd) \(t>i) = 9^171^, (0\ss\(f) s ) = g^m 2 ^, 
(<j> a \BB) = (2n) 4 5 4 (p^-p B +PB)9^BB, 

we find 

{MEr" = -(M^ « -2mg^ iBB - m a g^ mBB 

~ — (4m + m s )gg' + (m s — 2m) (B\Y\B) gg' , 



where in the second line we further assume the couplings to be given approximately by 
the additive quark model with SU(3) symmetry, so that g = g^ ~ g ( f >a while g^BB — 
(B\Y + 2\B)g' and g^ BB ~ (B\l - Y\B)g'. 

To estimate (M B ) opi we return to Weinberg's determination of the quark masses, 
which is performed assuming that the baryon matrix elements of ijjif) are proportional to 
the quark number, and is phenomenologically accurate to about 20%. Because all contri- 
butions to ipip such that this identification is correct come from one particle irreducible 
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diagrams [ 10 1 [ 11 ] , the successfulness of this approach may now be understood as a conse- 
quence of the cancellation eq. (|2.8|) , thus implying that (Mj^) opi is approximately flavor 



singletH. 

In view of this, we may conclude that the bulk of the flavor nonsinglet component of 
M% is provided by (M^) opr , and thus that Mf « (m a - 2m)gg' while Mf « 0. Thus 
we estimate AM in eqn. (|2.4|) to be very roughly of order 200 MeV if the couplings are 
of order unity, which is of the same order of magnitude, and the same sign, as the matrix 



element (N\Os\N) deduced from the mass splittings (|1.6|) . Consequently the theoretical 
estimate of as is approximately doubled, to around 50 MeV, which is quite compatible 
with the identification of as and a without the necessity for any violation of the Zweig 
rule. 

Whereas this estimate of the size of AM relies on the crude pole-dominance approx- 
imation ( |2.9| ) used to deduce eq. (|2.11|) , and on the neglect of (M|) opi , the failure of the 
identification of the mass splitting operator Os with Tis may be inferred from the exact 



results (|2.2|) , (|2.4| ) and (|2.8| ). Clearly all of the usual consequences of the assumption 
that baryon mass splittings are given by matrix elements of an octet operator Os, such 
as the Gell-Mann-Okubo mass formula, are intact; it is only the identification of Tis with 
that operator which is modified. The same applies to results derived assuming that an 
additive quark model picture applies to the matrix elements of ^0, such as Weinberg's 
determination of current quark masses || , since these only test the identification of baryon 
matrix elements of Os with (M^) opi which is correct due to the cancellation eq.( |2.8| ). The 
sigma term puzzle is thus resolved by the observation that the gluon condensate provides 



a positive contribution to the parameter M\ (which is negative) in eq. (|1.5|) ; the quark's 
contribution to the splittings is necessarily underestimated if this is not taken into account. 



3. Computation of Dynamical Quark Masses 

How can one test quantitatively the explanation of the sigma term puzzle proposed 
in the previous section? Because the gluon condensate provides the bulk of the baryons' 
masses, the SU(3)-dependence of its matrix elements correspond to a similar SU(3) de- 
pendence of the dynamical contribution to constituent quark masses. In particular, the 
dynamically generated mass should display a strong dependence on the current mass, an- 
ticorrelated to it. 



A discussion of the dynamical reason for this within the light-cone parton model is given in 
rei.pl. 
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3.1. The Quark Self Energy 



We can try to find the flavor dependence of the dynamically generated mass by study- 
ing the quark self energy, which we may compute by solving the quark Schwinger-Dyson 
equation [|T7f1 [|T9|| ; 



S-\p) = Z 2 S^{p) -Z l9 2 j ^A«V#(p - k )SmliP, k), (3.1) 

where S(p) = Z(p 2 )/(j/ + T>(p 2 )) is the full quark propagator, E(p 2 ) is the quark self 
energy, So(p) = (j/ + mo) -1 is the bare (or "current") quark propagator, D^Kp — k) the 
full gluon propagator, and T b v {j>, k) the full quark-gluon vertex function. Z\, Z 2 and Z 3 are 
the usual renormalization factors for the vertex, quark and gluon fields, respectively. The 
generators A a of the color group 577(3) are normalized as tr(A a A&) = -|<5 a b, and henceforth 
we will work throughout in Euclidean space. 

As it stands eq. (|3.1|) is not consistent with the renormalization group equations for 5, 
D and T, so it must be renormalization group improved ||20|| ; this results in the replacement 



of the coupling constant with a running coupling g 2 (p, k) inside the integral, where g(k,p) = 
g(p, k) and g{k,p) ~ g(k 2 ) when k 2 ^> p 2 , g(k 2 ) being the usual running coupling [|2T| [|i~9f . 

Let us now assume (neglecting ghost contributions) that r™ = A a r„, where Y v satisfies 
the Ward-Takahashi identity 

(p - k) v Y v (p, k) = S-\p) - S-\k). (3.2) 



Besides ensuring that Z\ = Z 2l eq.(|3.2j) may be used to find a consistent expression for 
the vertex in terms of the quark propagator. A suitable solution was given long ago by 



Landau [221: 



fc ) = fh^f ~ S ~ 1{ ^ + TflAP ~ k) ^ 2 ^ (3 - 3) 

where = (5^ — p^Pv/p 2 ), Z{p,k) = Z(k,p) and Z(p,k) ~ Z(k 2 ) for k 2 ^> p 2 . Even 
though this ansatz has a kinematic singularity in the infrared (i.e. as p — > k), it is sufficient 
for our purposes since it guarantees multiplicative renormalizability, and thus will give a 
quark self energy which has the correct gauge independent asymptotic behaviour at large 
p 2 (see ref. [20| for a more complete discussion). This confines ambiguities to the infrared 
region over which we have little control anyway. Analogously, the gluon propagator can 
be taken to have its asymptotic form, up to a momentum-dependent function d(p) which 
parametrizes infrared uncertainties; 



D?M = 5 ah d{p) T ^ v)+ f^ ,P \ (3.4) 



where d(p) — > 1 as p 2 — > oo. 
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Using the vertex ( pT3|) and the propagator fl3.4|), the Schwinger-Dyson equation ( |3TT| ) 
becomes 



S-\p) = Z 2 So 1 (p)+SC 2 Z 2 



d 4 k g 2 (p,k)d(p — k) 
(2<F (p-k) 2 



x 



T^(p - k) lfl S{k) lv Z{p, k)- 1 + ijt^. (S{k)S-\p) - 1) 



(3.5) 



where the colour factor C 2 



2N r . 
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Eq. (|3.5|) can be further simplified by assuming that the a priori unknown factors 
g(p,k), d(p — k), and Z(p,k) can each be approximated by the asymptotic values they 
take when one of their arguments is much larger then the other, i.e. g(p, k) 2 d(p — k) ^ 
g(max(p 2 , k 2 )) 2 and Z(p,k) ~ Z(max(p 2 , A: 2 )) [p2| . We may then perform the angular 
integrals analytically to give 



Z~\p 2 



z 2 + 



iz 2 c 2 

16tt 2 



g 2 (p 2 ) E(p 2 



p 2 k 2 z{k 2 mk 2 ) 

P" Z(p 2 ) Jo (k 2 + E 2 (k 2 )) 

A 2 ~2(U2\ 



dk 



2 g\k 2 )Z(k< 



X(p 2 ) = Z(p 2 )Z 2 m (A) + 



3C 2 
16tt 5 



£; 2 + E 2 (£; 2 )_ 
9 2 (P 2 _1 

Jo 



+ Z(p 2 



A 



2 <7 2 (fc 2 )E(fc 2 ) 
2 & 2 + E 2 (£; 2 ) 



" 2 fc 2 Z(fc 2 )E(fc 2 ) 
" £; 2 + E 2 (fc 2 ) 

2 



+ ^E(p 2 ) / rffc 

Jv 2 



2 g_^)Z{y 
k 2 + E 2 (k 2 



(3.6) 



where we have introduced an ultraviolet momentum cutoff A. In the Landau gauge £ = 0, 
Z(p 2 ) = Z 2 , and it is not difficult to see that ( j3lf ) reduces to the same equation as 
obtained in the (renormalization group improved) ladder approximation |19|1 . 

Eq. (|3.6|) is sufficiently simple that it may be solved numerically for the quark self 
energy E(p 2 ). All infrared uncertainties (in the regions k 2 ~ M 2 , p 2 ~ M 2 and (p — k) 2 ~ 
M 2 , where M 2 is the "QCD scale") have been rather crudely absorbed into the infrared 
uncertainty in the form of the running coupling g(p 2 ). By taking suitable parameterizations 
of this uncertainty we hope to get some feel for the overall uncertainty due to our ignorance 
of the physics of the infrared region. We may assume for example that 



1 



f3 lii{6 + p 2 /M 2 ) 



(3.7) 



where /3q = 11N / R 1 Nf , and 5 may be adjusted at will to give various values of the strong 



48tt 

coupling at p 2 = M 2 
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3.2. Current Mass and Dynamical Mass 

After performing such a calculation, it is still necessary to adopt some procedure to 
extract the quark "mass" from its self energy. Since we only have a rather approximate 
form for the self energy in Euclidean space, there is no question of being able to continue 
to Minkowski space to search for a pole in the quark propagator, even if one believed that 
such a pole should exist. We may however estimate the region of p 2 in which the self 
energy makes its dominant contribution to the mass of an on-shell baryon. The simplest 
such estimate is to take the "pole-mass" defined as the solution to ms = £(raf}): the self 
energy is taken to yield the mass when evaluated at the scale set by the mass itself. 

We may assess the uncertainty in this estimate by observing that the baryon mass 
could be determined, if we knew its Bethe-Salpeter amplitude, by calculating an integral 
over the relative momenta of its constituent quarks; the effect of these integrals would be 
to give the contribution of each quark as some weighted mean of its self energy. We may 
attempt to simulate this by smearing the quark's self-energy with a momentum-dependent 
"form factor" p(p 2 ), and then defining the mass as the solution of 

I* oo 

m£ = (E> p = / dp 2 p(p 2 )E(p 2 ), (3.8) 
Jo 

where p satisfies the normalization condition (l) p = 1 and the additional condition (p 2 ) p = 
(m^) 2 , which ensures that the pole-mass definition is reproduced when p tends to a Dirac 
delta. A suitable form of p could be p(p 2 ) = Afp 2 (p 2 + a) _n , for n = 4, 5, 6, say S ; with M 
and a determined by the two conditions. 

The masses determined in this way may be thought of as "constituent" masses; 
before we can use them to estimate AM we must still separate out the renormalized 
"current" mass mc from the dynamical mass which arises nonperturbatively to break 
chiral symmetry even when mc is small. This can be done by observing that the amount 
of explicit symmetry breaking may be extracted from the asymptotic behaviour at large 



p of the self energy [17 



S(p 2 ) ~ mln-\p 2 /M 2 )(l + 0(g(p 2 ) 2 ), (3.9) 



where A = 7^//3o = 2 n JtTn -2N f ) = I wneri N c = Nf = 3, and 7^ is the leading coeffi- 
cient in the perturbative expansion of the mass anomalous dimension. The renormalization 
group invariant mass parameter m is related to the running current mass in the usual way: 

p) ~ mln~ X (p 2 /M 2 ), (3.10) 

/j, — * 00 

where p is the renormalization scale. 

By varying the bare mass mo (A) we obtain thus a family of solutions to (|3.6|) , each 



of which corresponds to a different (unique) value of m. Calculating the constituent mass 



3 Both the quark propagators and the Bethe-Salpeter amplitude fall as powers when p 2 be- 
comes large. 
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raj: for each of these solutions gives us as a function of m. The asymptotic behaviour 
( |3.9|) in combination with the definition of tells us immediately that 



\rr A {m z /M I ){l + 0(g{m z Y)). (3.11) 

m — > oo 

To obtain the dynamical mass mii(m) we should now subtract from mj(m) the current 
mass mc eq.( [3.10|) , evaluated at the scale of the constituent mass: 

mofm) = ms(m) - mc(m, msfm)). (3-12) 



Due to the asymptotic behaviors of ms, eq.( |3.11| ), and mc, eq.( [3.10|) , the dynamical 



mass defined by eq. (|3.12|) vanishes for large m, as it ought to (the current mass and the 



constituent mass coincide for large values of m). 

It is still necessary however to further specify the form of mc in the infraredH Clearly 
ms(0) = mi)(0), so mc(0, mo(0)) = 0. Furthermore we may use Weinberg's definition || 
of the constituent quark mass, i.e., assume that in the chiral limit the dependence of mo 
on m is dominated by the linearized dependence of mc on m (which, recalling eq. ( |2.3|) , is 
essentially the dependence of Mg on the quark content of the given state). In other words 
we assume that for small m 

ms(m) - mfl(0) + mc( m ) weM), (3.13) 
and, identifying the linear terms in a Taylor expansion in m, 



dm 



c 



dm 



dm-y 



(3.14) 

m=0 



Note that this condition may be enforced only at one particular m since m^{m) is nonlinear. 
In practice we choose a smooth interpolating function to parameterize mc for all m, such 
as, for example (compare eq. (|3.7|) ) 

mc(m, ms(m)) = m ln _A (e + m 2 /M 2 ), (3.15) 



where the parameter e is chosen to ensure the satisfaction of ( |3.1 1 ) 



The dynamical mass m£>(m) is then a smooth function, flat at the origin, tending 
gradually to zero on a scale of mo{0) ~ 300 MeV. It is this mass which should be 
identified with the dynamical contribution (see (|2.3|) ) at the constituent quark level, 
since the "current" quark contribution is clearly already accounted for by mc- Thus, 
referring back to (|1.5[ ) and (|2.4| ) we see that assuming a linear mass formula 



AM ~ Mf ~ m D (0) -m D {m s ), (3.16) 

where m s is the value of m appropriate for the strange quark mass; taking the running 
mass m c (m s , 1 GeV) = 175 ± 50 MeV [§ we find m s = 200 ± 50 MeV [§, which is of 
the same order of magnitude as mo(0). From the qualitative behaviour of m^im) we 
thus expect AM to be loosely of order 150MeV, in rough agreement with the estimate of 
section 2.2. 



Even if we could compute mc (m, fi) reliably for small /x, we would still require extra condi- 
tions to fix the (probably rather large) scheme dependence. 
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3.3. Numerical Results 

We may now compute AM eq.( |3.16|) by solving the Schwinger-Dyson equation in 
the approximate form ( |3.6| ), with N c = Nf = 3 and, for example, the ansatz ( |3.7|) , and 
extracting m£>(ra) according to (|3.8|), (|3.12| ) and ( |3.15| ). The results of such calculations 



with a s (0) = (47r/?o ln5) _1 = 2.5 and the gauge parameter £ = are displayed in fig. 2, for 
both the "pole" definition and the various averaging definitions ( |3.8|) of the constituent 
mass; it can be seen that m_D(m) depends very little on which definition we adopt, so all 
remaining calculations are done with the simpler "pole" definition. 

The variation with gauge parameter £ is explored in fig. 3 — again we find that for 
|£| ~ 1 mo is approximately independent of the choice of gauge, as expected since the 
Landau vertex (|3.3| ) satisfies the Ward-Takahashi identity (|3.2j ) . Of course for larger values 
of £ there is still some gauge dependence, resulting from the infrared uncertainties. This 
could presumably only be tamed by using a vertex which satisfied the full Slavnov-Taylor 
identity, and thus including ghost contributions which are significant in the infrared, but 
this is not possible with the present level of expertise!. 

However the dominant uncertainty comes not from gauge dependence but from our 
general ignorance of the infrared. If we vary the strong coupling in the infrared by choosing 
different values for S in ( |3.7|) , we can obtain wide variations in m_o(m); displayed in fig. 4 
are curves for 5 = 2, 1.5, 1.2, 1.01 corresponding to a s (0) = 0.5, 1, 2.5, 50. Clearly as a s (0) 
is increased the dynamical mass becomes firmer, and consequently falls off less rapidly as 
m is increased; AM is thus reduced. To check this observation we also computed curves 
for couplings containing various infrared singularities, for example 5 A (p) or l/(p 2 ) a with 
< a < 1 (for a = 1 the Schwinger-Dyson equation is infrared divergent) which confirm 
this expectations: the stronger the singularity, the firmer the dynamical mass. 

We conclude that it is not possible to compute AM without further (highly nontrivial) 
information on the infrared structure of QCD — all we can do is estimate AM to lie in 
the range 50 MeV ^ AM 250MeV which translates (using (|2.4[) ) into the estimate 

30 MeV ^ a 8 ^ 60 MeV. (3.17) 

in perfectly satisfactory agreement with the value of a inferred from experiment (and with 
the estimate presented in section 2.2). 

Amusingly, we also have sufficient information to compute a using the so-called 
Feynman-Hellman formula at the constituent quark level, which amounts to the (rather 
reasonable) assumption that mc is approximately linear in the small light flavor average 
mass rh: 

a ~ 3m— — = 3mc(fn, mnifi)) (3.18) 
dm 



where in the last step we used eqns. (|3TT3|) , (|3.14|) , and we may compute mc(m,m£i(0)) 



through eq.( |3.15| ) from the known value of the running current mass at the nucleon scale. 
Taking m c (fh, 1 GeV) = 5 ± 3 MeV [§,0, we find a 8 = 40 ± 20 ± 15MeV, the first error 



It is however possible to do computations with more complicated vertices |23| which satisfy 



the Ward-Takahashi identity (3.2) but are however free of the infrared kinematic singularity which 



plagues fl3.3|) . The results are however little different to those with the simpler Landau vertex. 
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being due to the uncertainty in mc(m, 1 GeV), and the second due to the uncertainty in 
e, which we estimate by varying a s (0) from 1 to 50. Since the former is rather large, and 
the latter probably underestimated, this is clearly a less reliable result than ( |3.17|) , but 
nonetheless it provides a useful independent check on the consistency of the calculation. 

The calculations presented here are necessarily extremely uncertain, not only due to 
our implicit use of a constituent quark picture of the baryon (instead of a linear mass 
formula, we could have used a more sophisticated one (see for example in which 

MrP would also be nonzero), but also to uncertainties in the separation of the dynamical 
contribution from the "current" quark contribution, in the form of the quark-gluon vertex 
function and, most importantly, our complete ignorance of the behaviour of the gluon 
propagator in the infrared. Note that the possible effects of quark loops, both light and 
heavy, are only a part of this uncertainty, in that while they certainly contribute to the 
gluon propagator in the infrared, there is no reason to believe that they dominate it. 

A more accurate computation of the sigma term than that presented here would 
thus require much more powerful techniques (for example lattice studies), because of its 
sensitivity to the infrared dynamics of chiral symmetry breakdown. Certainly nothing 
more can be learnt from models such as that of Nambu-Jona-Lasinio, or the chiral quark 
model, since these necessarily trivialize this dynamical structure by assuming a constant 
quark self energy. 

4. Conclusion 

In this paper we have shown that the traditional argument relating the sigma term 
to baryon mass splittings is incorrect, due to the neglect of the dynamical contribution 
to such mass splittings from the trace anomaly. This contribution may be viewed as a 
nonsinglet gluonic contribution to constituent masses, and as such it does not spoil all the 
usual quark model results, although it does contradict the naive assumption that mass 
splittings at the constituent and current level are the same. 

We have further attempted a computation of this contribution through the solution of 
the Schwinger-Dyson equation for the quark self energy, and find that it not only has the 
right sign, but is also quite large enough to agree with the experimental determination. 
However since it is rather sensitive to the infrared dynamics, it seems to us to be rather 
difficult to obtain a more precise determination, and thus a firm bound on the amount 
of Zweig rule violation, without a better understanding of the details of the mechanism 
responsible for the dynamical breaking of chiral symmetry. What we can say, however, is 
that without such an understanding the experimental value of the sigma term cannot be 
used as evidence of strange quarks in the nucleon. 

Rather than providing evidence for the breakdown of the naive quark model, the 
sigma term appears thus to provide nontrivial information on the dynamics of QCD in the 
infrared. In view of this, both a better experimental determination of the sigma term, and 
a better theoretical understanding of nonsinglet contributions to quark masses would be 
highly desirable. 
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Figure Captions 



1] The Schwinger-Dyson equation. 

2] The dynamical mass m£>(m) vs. the current mass m, both in units of mo(0) (solid); 
also m^(m) (defined in eqn.(|3.8|)) for n=4,5,6 (long dashes, short dashes, dots respec- 
tively) . 

3] As fig. 2, but for the gauge parameter £ = 0, 1, 2, — 1, —2 (solid, long dashes, short 
dashes, dots, dots and dashes respectively). 

4] As fig. 2, but for a s (0) = 0.5, 1, 2.5,50 (dots, short dashes, solid, long dashes respec- 
tively) . 
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